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Abstract 

An often used mass formula for Nambu-Goldstone (NG) bosons in QCD, 
such as the pions, involves the condensate < qq >, f n and the quark 
current masses. We argue, within the context of the Global Colour 
Model to QCD, that this expression is wrong. Analysis of the interplay 
between the Dyson-Schwinger equation for the constituent quark effect 
and the Bethe-Salpeter equation for the NG boson results in a new mass 
formula. 
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The computation of the low energy properties of QCD is a difficult non-perturbative 
problem in quantum field theory. However one species of hadron, the (almost) Nambu- 
Goldstone (NG) bosons, such as the pions, have always played a key role. Because they 
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are directly associated with the dynamical breaking of chiral symmetry their properties 
are strong indicators of the nature of the underlying quark-gluon dynamics in QCD. 
However the small current masses of the u and d quarks means that the pions are not 
strictly massless, as for true NG bosons, but acquire small masses. They are therefore 
also significant as low mass hadronic excitations. 

One expects that there should be some perturbative expression for the pion mass in 
terms of the quark current masses which is built upon the underlying non-perturbative 
chiral-limit quark-gluon dynamics. While the relation of the low pion mass to the 
breaking of chiral symmetry dates back to the current algebra era and PCAC [ij, the 
often used implementation in QCD has the form, 

Ml = (m “ + mj> " ( 1 ) 

J 7T 

where the integral p =<qq > is the so called condensate parameter. For N c = 3 

p = N c tr(G(x = 0)) = 12 J ^-^(g 2 ), (2) 

and f n is the usual pion decay constant. In (2) <x s (s) is the chiral limit ( m —> 0) scalar 
part of the constituent quark propagator 

G(q) = (LA(s; m)q.j + B(s ; m) + m)^ 1 = —iq.^a v (s\ m) + <r s (s; m). (3) 

We note that the expression for p in (2) is divergent in QCD, because for large 
s —> oo B(s) decreases like l/sln[s/A 2 ] 1-A where A = 12/(33 — 2 Nf) and A is the QCD 
scale parameter. Some integration cutoff is usually introduced. The values of m and 
< qq > are then usually quoted as being relative to some cutoff momentum, often 
1 GeV. An alternative approach Q is to use finite energy sum rules and Laplace sum 
rules. 

Here we present a new analysis of the chiral symmetry breaking in the Global 
Colour Model (GCM) to QCD. We extract a new expression for the pion mass, which 
essentially replaces (2). Nevertheless the new result, see equation (22) below, is very 
similar to (1) except that it contains a naturally arising cutoff function c(s), and also 
a dynamical enhancement function e s (s) for the quark current mass m. There is also 
a factor of 2 difference. In the Appendix we give an example of the type of incorrect 
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analysis that leads to the form in (2). The new mass formula implies that the pion 
mass is dominated by IR processes, and not UV processes as in (2). 

An overview and an insight into the nature of the non-perturbative low energy 
hadronic regime of QCD is provided by the functional integral hadronization of QCD 
§ |. This amounts to a dynamically determined change of functional integration 
variables, from quarks and gluons, to bare hadrons 

/ VqVqVAexp(-S qcd [A, q, q]+rjq + qrj) « 

j VTTVNVN...exp(-S had [Tr ,..., N, N,..} + J n [rj , rj\-K + ..) (4) 

The final functional integration over the hadrons gives the hadronic observables, 
and amounts to dressing each hadron by, mainly, lighter mesons. This functional in¬ 
tegral transformation cannot yet be done exactly. The basic insight is that the quark- 
gluon dynamics, on the LHS of (4), is fluctuation dominated, whereas the RHS is not, 
and for example the meson dressing of bare hadrons is known to be almost pertur¬ 
bative. In performing the change of variables essentially normal mode techniques are 
used 0. In practice this requires detailed numerical computation of the gluon propa¬ 
gator, quark propagators, and meson and baryon propagators. The mass-shell states 
of the latter are determined by covariant Bethe-Salpeter and Faddeev equations. The 
Faddeev computations are made feasible by using the diquark correlation propagators, 
which must also be determined. 

The first and easiest formal transformation results from doing the gluon integra¬ 
tions, leaving an action for quarks of the form 

S[q,q] = j q{x)(- 1 .d+M)q{x) + \ J ;j° J (x)j 0 Jy)D llu (x-y) + ^ J + - (5) 

where j“(x) = q{x)^- r y tl q{x). The GCM is a model field theory for QCD based on a 
truncation of S[q, q] in which the higher order n-point (n > 3) functions are neglected, 
and only the gluon 2-point function D^ u (x — y ) is retained. 

The GCM is thus a quantum field theory that can also be considered to be defined 
by the action 

Sgcm[q, q, = J ( q(x)(-'y.d + M + iA a J^-^^)5{x - y)q(y)+ 

+\ A li. x ) D ^(id)5{x - y)Al{y) S j (6) 
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where the matrix D~^ (p) is the inverse of D ( p ), which in turn is the Fourier transform 
of D^ u ^x). This action has a global colour symmetry. The GCM is thus analogous 
to QED except for colour currents and the non-quadratic phenomenological form for 
D~j}(p) in the pure gluon sector. The determination of the best-fit phenomenological 
form for D fiu (p) was reported in @ using the separable expansion technique to facilitate 
numerical computations and the fitting to some meson data. 

Having made the GCM truncation in (5) it is possible to proceed further and to 
transform Q the quark functional integrations into the hadronic functional integrations, 
as in (4). If the additional approximation D^i^x—y) —* g5 flL ,5(x—y) is made in (6), i.e. a 
contact coupling of the quark currents, then the NJL type models are obtained. If in (4) 
a derivative expansion of the complete non-local hadronic effective action is performed, 
then the Chiral Perturbation Theory (CPT) phenomenology is obtained. However 
in the GCM, with appropriate D^ u {x), all computations are finite and no cutoffs or 
renormalisation procedures are used. As well, using a mean field approximation, the 
soliton phenomenology for the baryons may be derived Q, and has been studied in [|6|]. 

In analogy with QED the GCM involves the determination of various coupled equa¬ 
tions for the propagators. In the approximation often employed (which is motivated 
by the functional integration analysis of (6), see j|]) the first equation is the Dyson- 
Schwinger (DSE) equation for the constituent quark propagator (the so-called vacuum 
equation of the GCM ||, |U ), 


B(p 2 -m) = 


16 


d 4 q 


D(p-q). 


B(q 2 ; m) + m 


q 2 A(q 2 -,m ) 2 + {B{q 2 \m) + m) 2 ’ 


(7) 


[A(p 2 ; m ) — 1 ]p 2 = 


d A q 


q.pD(p - q).- 


A(q 2 ;m) 


( 8 ) 


(27t) 4 1 y u q 2 A(q 2 - 1 m) 2 + (B(q 2 ] m) + m) 2 ’ 
using a Feynman-like gauge, D^ v (p) = 6 flu D(p), and the perturbative quark-gluon 
vertex function. The Landau gauge can also be used; see [|?| for comparison. 

Using Fourier transforms (7) may be written in the form, here for m = 0, 

3 B[x) 


D(x) = 


(9) 


16 a s (x) ’ 

which implies that knowledge of the quark propagator determines the effective GCM 
gluon propagator. Multiplying (9) by B(x)/D(x), and using Parseval’s identity for the 
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RHS, we obtain the identity 


/ 


d‘x Bix) 


D{x) 


16 

y 


A 


B(q)a s (q). 


( 10 ) 


The second basic equation is the Bethe-Salpeter equation (BSE) for the pion mass- 
shell state at the level of approximation that matches (7) and (8) in the GCM analysis 

0 , 1 

T f (p, P) = lf D (P ~ q)trsF(G + TVG-Tf)T3(q, P ) (11) 

where G± = G(±q — -^). This BSE is for isovector NG bosons, and only the dominant 
T = r^r -^75 amplitude is retained (see Q for discussion); the spin trace arises from 
projecting onto this dominant amplitude. Here {T b , b = 1 ,Np — 1} are the generators 
of SU(N f ), with tr(TfT 9 ) = \5 fg . 

The BSE (11) is an implicit equation for the mass shell P 2 = —M 2 . It has solutions 
only in the time-like region P 2 < 0. Fundamentally this is ensured by (7) and (8) being 
the specification of an absolute minima of an effective action after a bosonisation ||]. 
Nevertheless the loop momentum is kept in the space-like region q 2 > 0; this mixed 
metric device ensures that the quark and gluon propagators remain close to the real 
space-like region where they have been most thoroughly studied. Very little is known 
about these propagators in the time-like region q 2 < 0. The GCM gives a detailed 
description of the pion properties, including its coupling to other states, in the language 
of effective non-local actions 0]. 

The non-perturbative quark-gluon dynamics is expressed here in (7) and (8). Even 
when m = 0 eqn.(7) can have non-perturbative solutions with B ^ 0. This is the 
dynamical breaking of chiral symmetry. 

When m = 0 eqn.(ll) has a solution for P 2 = 0; the Goldstone theorem effect. For 
the zero linear momentum state {To = 0, P = 0} it is easily seen that eqn.(ll) reduces 
to eqn.(7) with T^(g, 0) = B(q 2 ). When P ^ 0 then T^(q,P) / B(q), and (11) must 
be solved for T^(q,P). 

We shall now determine an accurate expression for the mass of the pion when m ^ 
0. This amounts to finding an analytic solution to the BSE (11), when the constituent 
quark propagators are determined by (7) and (8). The result will be accurate to order 

m. 
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For small m ^ 0 we can introduce the Taylor expansions in m 


B(s ; m) + m = B(s ) + m.e s (s) + 0(m 2 ), (12) 

A(s, m) = A(s) + m.e v (s ) + 0{m 2 ). (13) 

For large space-like s we find that e s —> 1, but for small s we find that e s (s) can be 
significantly larger than 1 (see Fig.l). This is a dynamical enhancement of the quark 
current mass by gluon dressing in the infrared region. Even in the chiral limit the quark 
running mass M(s ) = B(s)/A(s) is essential for understanding any non-perturbative 
QCD quark effects. At s = 0.3GeF 2 we find @ that M(s) « 270 MeV. 

Because the pion mass M n is small when m is small, we can perform an expansion 
of the Pf, dependence in the kernel of (11). Since the analysis is Lorentz covariant we 
can, without loss of validity, choose to work in the rest frame with P = (■ iM n , 0), giving, 
for equal mass quarks for simplicity, 

r (p) = J D(p-q)I(s)T(q) + 

+ Y / j2^y D ( P ~ q) s(A(s) + e v (s).m) 2 + ( B(s ) + m.e s (s)y T<<q) + (U) 

where 

7(s) = 6 (cr 2 - 2{a s a' s + sa v a' v ) - s{a s a" - (a' s ) 2 ) - s 2 (a v a" - {a' v ) 2 ^j . (15) 


By using Fourier transforms the integral equation (14), now with explicit depen¬ 
dence on can be expressed in the form of a variational mass functional, 

M fri 2 =_— /' _IM!_ + 9 f 

7r[J U[TYJ {2^Ys(A{s) + e v {s).mf + (B{s)+m.e{s))^2U[TfJ D(x) 

(16) 

in which 

Mr ] 2 = f (04^)r(g) 2 . (it) 

The functional derivative 5M n [T] 2 /5T(q) = 0 reproduces (14). The mass functional 
(16) and its minimisation is equivalent to the pion BSE in the near chiral limit. To 
find an estimate for the minimum we need only note that the change in M 2 from its 
chiral limit value of zero will be of 1st order in m, while the change in the zero linear 
momentum frame T(q) from its chiral limit value B(q 2 ) will be of 2nd order in m. 
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Hence to lowest order in m we have that the pion mass is given by 

, , 2 _ 48m f d 4 q e s (s)B(s) + se v (s)A(s) B(s) 2 
~~ U[B] 2 J (2vr) 4 sA(.s) 2 + B{s) 2 sA(s) 2 + B(s ) 2 

24 r d A q B{s) 2 9 f 4 B(x) 2 2 

/,[5] 2 J (2vr) 4 sA (. S ) 2 + B(.s) 2 + 2/^R] 2 J ' X D(x) + ^ 

However the pion mass has been shown to be zero in the chiral limit. This is confirmed 
as the two O(m 0 ) terms in (18) cancel because of the identity (10). Note that it 
might appear that f n would contribute an extra m dependence from its kernel in (15). 
However because the numerator in (16) is already of order m, this extra contribution 
must be of higher order in m. 

Hence we finally arrive at the analytic expression, to 0(m), for the NG boson 
(mass) 2 from the solution of the BSE in (11) which includes the non-perturbative 
gluon dressing to give constituent quarks 



2 = 2mp eff 

TV f2 

J 7r 


(19) 


where 


f d A q 

Peff = 24 y 7^34 (e s (s)cr s (s) + se„(s)c7„(s)) c(s) 


( 20 ) 


defines an effective condensate parameter p e ff- The p e ff integrand involves a naturally 
arising function 


c(s) = 


B(sY 


sAisy + Bis) 2 ’ 


( 21 ) 


which acts as a smooth cutoff function. It is this function which causes the pion mass to 
be IR dominated. Of course the determination of B(s) and H(s) by (7) and (8) involves 
both IR and UV momenta. As well p e ff contains a contribution from the vector part of 
the chiral-limit quark propagator. Because this vector part makes a small contribution 
we can approximate (20) by the form 


Peff = 



d 4 q 

(2vr) 4 


e s (s)c(s)a s (s) 


( 22 ) 


This result contradicts (2). Significantly the factor c(s) makes (22) convergent when 
B(s ) has the QCD determined asymptotic form discussed after (3). The e enhance¬ 
ment function is important to the numerical values but does not affect the convergence 
analysis because asymptotically lim is ^ 00 e s (s) = 1. 
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We now give examples of the functions occuring in (22). In ref. j?]] the GCM phe¬ 
nomenological chiral-limit quark propagator and the corresponding gluon propagator 
were obtained by fitting a number of meson observables to meson data. From these 
forms the cutoff function c(s) and the enhancement function e s (s) can be computed. 
The resulting functions are shown in Fig.l. We note that the value of e s (s ) ~ 4 implies 
that in the infrared region, appropriate to the internal dynamics of hadrons, the quark 
current mass of ~ QMeV is enhanced by gluon dressing to some 24MeU. Of course the 
major effect is the chiral-limit constituent mass of some 270MeV. We compute using 
(22) that p e ff = (0.233GeU) 3 , and from (17) that f n = 93.0 MeV. With m = 6.5 MeV 
(22) gives Mjf = 138.5IWeU without cutoffs or renornralisation procedures. We note 
that these numerical results are actually close to the usual values using the forms in (1) 
and (2) with the traditional cutoff of 1 GeV. This may be incidental to the particular 
parameterisation used for the phenomenologcal GCM propagators used in 0. 

It would be interesting to see if the corrected form in (22) can be rigorously obtained 
within QCD. 

Appendix: Here we present one form of the flawed analysis that leads to the incorrect 
result in (2). Consider the mass functional in (16). Suppose we make the ansatz that 

T n (q) = B(q 2 ;m). (23) 


From (7) we obtain by Fourier transforms and Parseval’s identity, in analogy to (10), 


’ 4 B(x-m) 2 

“dFT 


16 

T 


d A q 

w? 


B(q;m) 


B(q ; m) + m 


q 2 A(q\ m ) 2 + ( B(q ; m) + m ) 2 


Using (23) and (24) the mass functional (16) then gives 


2 24 

Ml = T7T 


d A q 


B(q\m) 


J (27r) 4 s(A(s) + e v (s).m) 2 + (B(s) + m.e(s)) 2 

24 r d A q B(q] m)(B(q; m) + m) 
f 2 J (27r) 4 q 2 A(q-, m) 2 + (B(q; m) + m) 2 ' 

24 I" d A q B(q\m)m 

f 2 J (27t) 4 q 2 A(q; m) 2 + (B(q; m) + m) 2 

24 f d A q B(q)m 


+ 


(24) 


(25) 

(26) 

“ f 2 J (27r) 4 q 2 A(q) 2 + (B(q)) 2 (27) 

to 0(m), and this gives eqns(l) and (2). The flaw in this argument is to overlook that 
in perturbation theory eigenvalues change at 1st order and eigenvectors change at 2nd 









order of the perturbation parameter. The ansatz (23) actually gives T 7r ((/) a hidden 
1st order dependence on m, and is thus incorrect. Essentially the same error has been 
made in the functional integral formulation of the GCM, where the ansatz (23) arises 
through the choice of the non-local meson fluctuation. In all GCM papers p should be 
replaced by p e ff- 

RTC acknowledges useful conversations with C. J. Burden, C. D. Roberts, P. C. 
Tandy and M. Thomson. C. D. Roberts [|9] has independently arrived at similar results 
to those reported here. 


Figure Caption 

Fig.l Shows (a) the scalar enhancement function e s (s), (b) the integrand se s (s)c(s)cr s (s) 
(from eqn.(22)) of the effective condensate parameter p e ff (in arbitrary units), and (c) 
the cutoff function c(s). 
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